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ABSTRACT 


This  report  summarizes  the  activities  carried  out  under  Signal 
Corps  Contract  No.  DAAB07-68-C-0381  (Texas  A&M  Research  Foundation 
Project  586)  during  the  contract  period  15  June  1969  through 
15  December  1969. 

Activities  during  this  period  have  centered  on  preparation  of 
the  equations  for  the  forested  boundary  layer  for  solution 'on 
the  general  purpose  analog  computer  at  Texas  A&M  University. 

These  expressions,  which  were  presented  in  Semi-Annual  Report  2, 
have  been  programmed  for  a  seven-level  simulation,  three  in  the 
forest  and  four  in  the  free  air  above,  and  for  the  surface  and 
free  air-canopy  interfaces.  The  wiring  of  the  patchboards  is 
proceeding  as  the  wiring  diagrams  are  completed. 

Some  qualitative  diffusion  tests  were  carried  out  in  the  tropical 
forest  in  the  Panama  Canal  Zone.  Visual  and  photographic  observa¬ 
tions  have  been  summarized. 
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I.  Introduction 

An  experimental  system  of  equations  has  been  developed  to 
simulate  the  atmosphere  within  a  forest  and  in  the  atmospheric 
boundary  layer  just  above  the  forest.  These  equations,  accom¬ 
panied  by  a  discussion  of  the  physical  considerations  used  in 
the  development  of  the  set,  were  Riven  in  Semi-Annual  Report  2. 

They  are  repented  in  Appendix  A  for  convenient  referral.  Except 
for  some  minor  changes  in  symbology,  the  only  modification  has  been 
in  the  calculation  of  stress  at  the  top  of  the  canopy.  In  the  new 
set,  a  logarithmic  distribution  of  wind  speed,  rather  than  a 
linear  one,  has  been  assumed. 

The  essential  parts  of  the  formulation  are  the  equations 
which  define  the  temporal  variations  of  momentum,  temperature, 
and  vapor  pressure  in  the  forest  and  in  the  free  air  above  (equations 
lf-4f,  1-4)  and  those  which  represent  the  heat  balance  at  the  top 
of  the  canopy  and  at  the  earth’s  surface  (equations  IQ,  20,  1.8f,  l«f) 
Host  of  the  other  expressions  define  the  various  terms  in  these 
basic  equations  or  are  approximations  which  simplify  their  solution. 
Many  of  these  are  standard,  e.g.,  the  relationship  between  specific 
humidity  and  vapor  pressure  (equation  11) ;  others  are  specific  for 
this  particular  problem  as  in  the  expression  for  the  momentum 
exchange  coefficient  in  the  free  air-canopy  section  (equation  37). 

The  principle  activities  over  the  past  six  months  have  been 
devoted  to  the  task  of  programming  the  eouations  for  solution 
on  the  general  purpose  analog  computer  at  Texas  A&M.  The 
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scaled  analog  format  which  is  being  wired  for  the  simulation  is 
discussed  in  section  II  of  this  report.  Section  III  contains 
a  summary  of  the  visual  and  photographic  observations  made  during 
some  oualitative  diffusion  tests  carried  out  in  the  Panama  Canal 


Zone  forest  during  August. 


II.  Scaled  Analog  Format 

The  general  analog  form  for  those  equations  which  will  be 
solved  on  the  computer  are  given  in  Appendix  B.  Relationships 
such  as  equations  (11)  and  (13)  which  define  terms  in  other,  more 
basic,  equations  or  which  are  used  to  simplify  them  are  not 
included  in  the  listing.  The  numbering  used  in  Appendix  A  is 
retained  in  Appendix  B  with  the  prefix  "a1'  to  indicate  analog 
format.  Thus,  equation  (a7)  is  the  analog  equivalent  of  the 
mathematical  form,  equation  (7),  in  Appendix  A. 

The  transformation  from  mathematical  to  analog  format  is 
effected  primarily  by  rhe  integration  of  the  equation.  Where 
necessarv,  simpler  or  more  readily  solved  forms  of  the  individual 
terms  are  substituted  for  the  original  ones.  Thus  in  equation  (al) 
the  value  of  the  x-componer.t  of  the  wind  at  any  time,  t,  is  given 
by  the  integration,  from,  initial  time  to  time  t,  of  the  several 
terns  on  the  right.  The  initial  conditions  have  been  omitted  as 
understood,  but,  of  course,  will  be  included  in  the  simulation. 
Finite  differences  have  been  used  as  an  approximation  of  differ¬ 
entiation  in  the  vertical,  usually  with  the  implicit  assumption 
of  linear  variation  of  the  meteorological  parameter  between  the 
grid  points.  The  individual  terms  in  each  equation  will  be 
discussed  below,  after  the  scaled  equations  are  given. 

Subscript  indices  are  used  to  indicate  the  height  at  which  the 
various  parameters  are  computed.  The  subscript  j  signifies  the 
level  at  which  the  simulation  of  the  primary  variables  is  carried 
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out.  The  bottom  of  the  1-yer  containing  the  j**1  level  is 
signifiel  by  the  double  subscript  "ij"  and  the  top  by  the  double 
subscript  "jk".  Thus,  the  jth  level,  at  height  x^,  lies  between 
height  z^  at  the  bottom  of  the  layer  and  height  z^  at  the  top 
of  the  jtl(  layer.  The  depth  of  that  layer  is  indicated  by  Az^. 

The  bottom  of  the  layer  for  the  j**1  level  coincides  with  the  top 
of  the  layer  for  the  i  level  (i«j-l).  Thus,  x  . .  is  the  x- 

3C,lj 

component  of  the  shearing  stress  at  the  interface  of  the  layers 
containing  the  i^  and  levels.  The  reader  is  referred  to 
figure  1  for  further  clarification  of  the  concepts  of  levels  and 
layers  and  the  usage  of  the  indices  "j"  and  "ij".  The  primary 
variables,  u,  v,  T,  and  e,  are  predicted  for  the  levels  j  while 
the  stresses  and  fluxes  are  computed  at  the  interface  c  ?  the  layers 
containing  successive  levels,  that  is,  at  heights  Indicated  by 
the  double  index  ij. 

In  preparing  for  the  simulation  each  of  the  equations  in 
Appendix  B  was  expanded  to  account  for  the  multiplicity  of  levels. 

The  expanded  set  has  been  scaled  so  that  the  voltage  representing 
each  variable  will  cover  as  large  a  range  as  feasible  during  the 
simulation  without  overloading  the  amplifiers.  The  standard 
scaling  for  the  general  purpose  analog  computer  at  Texas  A&M  is 
based  on  *  reference  unit  of  100  volts.  The  scale  factors  which 
have  been  adopted  are  estimates  of  the  maximum  values  (absolute) 
which  the  meteorological  variables  are  likely  to  attain  in  this 
simulation.  These  are  represented  on  the  computer  as  100  volt  outputs. 
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Figure  1.  Schematic  presentation  of  the  levels  and 
layers  of  the  forest  model  to  be  used  for  the  initial  tests. 
Levels  are  indicated  by  solid  lines  and  layer  interfaces  by 
dashed  lines.  The  numbers  in  the  center  *re  the  indices  used 
to  signify  the  levels  and  layer  Interfaces  at  heights  given 
along  the  edges.  The  top  of  the  canopy  is  at  40  m. 


Tentative  sraling  factors  for  all  variables  are  listed  in 
Table  1,  along  with  certain  other  machine  constants.  These  are 
subject  to  change  should  the  test  solutions  indicate  serious 
over-  or  under-scaling  of  the  parameters.  The  tine  scaling,  computer 
time  versus  real  time,  of  1440  to  1  will  be  used  in  the  initial 
simulation.  That  is,  the  computer  solution  will  he  1440  times 
faster  than  the  real  problem  time,  so  that  an  eight-hour  predictions 
can  he  obtained  in  about  20  seconds. 

The  scaled  equations  are  given  in  Table  2.  They  are  numbered 
to  agree  with  the  corresponding  mathematical  form  in  Appendix  A, 
with  the  prefix  'V  to  indicate  scaled  form.  Thus,  the  scaled 
form  of  the  mathematical  expression  for  the  x-component  of  the 
shearing  stress  (equation  7,  Appendix  A)  is  equation  (s7)  in 
Table  2.  just  as  the  analog  form  is  equation  (a7)  in  Appendix  B. 

In  the  scaled,  representation  the  forested  boundary  layer  is 
treated  in  seven  layers,  three  in  the  forest  itself  and  four  in 
the  free  air  above  the  canopy  (figure  1).  The  wind  components, 
temperature  and  vapor  pressure  at  the  four  levels  above  the 
forest,  are  found  from  equations  (si)  -  (s4)  and  those  at  the 
three  levels  in  the  forest  in  equations  (slf)  -  (s4f).  The  stresses 
and  the  fluxes  of  heat  and  momentum  at  the  layer  Interfaces  are 
given  in  equations  (s7)  -  (slO)  for  the  region  above  the  canopy  and 
in  equations  (s7f),  (s8f),  (sllf),  and  (sl2f)  for  the  forest 
atmosphere.  In  addition  the  equations  which  provide  the  temperatures 
and  vapor  pressures,  the  stresses,  end  the  heat  and  vapor  fluxes, 
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Table  1.  Scale  factor  and  constants  employed  in  the  forest 
boundary  layer  simulation. 


A.  Scale  Factors 


e 

F*»Fy 

FCG 

ICG 

K, 

q0 

q« 


10  cm/sec 
100  mb 
1000  cm/ sec3 
100  m/100  km 
100  m/100  km 
50,000  cn?/sec 
0.05  cal/cnPsec 
0.05  cal/cmPsec 


Rn  :  0.05  cal/cnfsec 

S^  :  20  cm/sec 

S,/  :  10  cra/sec 

A 

T  :  50 °C 

time  :  1440:1 

u,v  :  50  m/sec  (free  air) 

:  10  m/sec  (forest) 

T*,Ty:  20  dynes /cn? 


B.  Constants 


0.6210 

P 

»  86,400  sec 

0.2400  cal/gm  deg 

Po 

-  1012  mb 

980  cm/ sec3 

Z 

-  1050  m 

0.0330  cal/co^sec 

CT 

•  1.354xl0"13  cal/caPsec  deg 

0.4000 

U) 

■  7 . 3xl0~5  rad/ sec 

593.3  cal/gm 

P 

»  1 . 2xl0-3  gm/cm9 
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Table  2.  Scaled 
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at  the  surface  and  at  the  top  of  the  canopy,  are  given  for  these 
very  special  regions  in  the  sections  indicated  by  B  and  D,  (equations 
sl9  -  s37  and  sl8f  -  s37f  respectively). 

The  specific  levels  and  layer  depths  used  in  a  simulation 
do  not  enter  into  the  scaled  equations  except  in  the  coefficients 
which  are  combinations  of  geometric  dimensions,  non-time-dependent 
physical  variables  and  factors  required  to  maintain  algebraic 
correctness  when  scaling  the  meteorological  parameters.  The 
depth  of  forcat ,  d,  to  bo  used  in  the  lntiol  tests  is  40  meters 
and  the  depths  of  the  surface  layer  (A1)  and  the  air-canopy 
layer  (2A)  are  1.5  and  10  meters,  respectively.  The  mathematical 
forms  of  the  coefficients  and  the  numerical  values  to  be  used  in 
the  initial  tent  solutions  arc  given  in  Table  3. 

The  wind  components  have  been  scaled  for  maxima  of  50  m/sec 
above  the  canopy.  Obser  ;at *.ns  Indicate  that  the  winds  in  the 
forest  are  quite  light  and  50  m/sec  would  probably  be  a  gross 
over-estimate  of  the  maximum  values  likely  to  be  attained. 

Therefore,  the  wind  components  in  the  forest  have  been  scaled 
for  10  m/sec.  The  temperature  is  scaled  for  50°C  and  the  vapor 
pressure  for  100  mb  both  above  and  within  the  forest.  Thus,  the 
principal  computer  variables  for  which  solutions  are  sought 
become  u/5000,  v/5000  in  the  free  air,  u/1000,  v/1000  in  the 

forest,  T/50  and  e/100  both  in  the  free  air  and  in  the  forest. 

2 

The  eddy-stresses  have  been  scaled  to  20  dynes  per  cm  ,  the 

2 

sensible  and  latent  heat  fluxe6  to  .05  cal/cm  sec.  The  scale 
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Table  3.  Form  and  magnitude  of  the  coefficients  in  the  scaled  format  of  the  equations  for  the  forested  boundary 
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Coefficient 


factors  employed  for  the  other  variables  will  he  covered  in  the 
course  of  the  discussion  of  the  individual  terms  of  the  scaled 
equations. 

A  one-to-one  correspondence  exists  between  the  general  and 
scaled  formats  of  the  same  number.  However,  it  is  not  always 
clear  how  the  one  follows  from  the  other,  partly  because  of 
the  level  and  layer  multiplicity  but  largely  because  of  the 
substitutions  from  the  defining  relationships  and  certain  approxima 
tions  which  simplify  the  solution.  Conseciuentlv ,  each  of  the 
terms  in  the  equations  for  the  basic  meteorological  parameters  are 
discussed  in  some  detail  below. 

A.  The  Wind  Components  (equations  si,  s2,  slf,  s2f) 

The  forms  of  the  equations  for  the  free  air  and  the  forest 
atmospheres  are  essentially  the  same  except  for  the  addition  of 
a  drag  force  in  those  for  the  forest.  It  is  fairly  simple 
to  follow  the  transformation  from  the  defining  mathematical 
expression  to  the  general  analog  format  except  for  the  advection 
terms.  These  have  been  re-expressed,  partially,  in  terms  of 
derivatives  obtainable  from  kinematic  analyses.  The  factor  of 
5  appearing  with  some  of  the  terms  in  the  forest  equations 
(slf,  s2f)  enters  because  the  coefficient  corrects  for  a  wind 
scale  factor  of  5000  cm/sec  whereas  the  winds  in  the  forest  are 
scaled  for  1000  cm/sec.  A  term-by-term  discussion  of  the  scaled 
east-west  wind  expressions  follows. 
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1.  The  Contour  Gradient 

The  contour  gradient  forces,  g  g  are  effectively 
scaled  by  the  time  and  velocity  scale  factors  and  the  need 
to  maintain  algebraic  correctness.  The  rather  peculiar 
scale  factor  of  3.473  cn/sec^  is  the  ratio  of  the  velocity 
and  time  scales  in  the  free  air. 

The  contour  gradient  at  any  level  is  obtained  through  a 
"buildup"  from  the  surface  value  in  successive  steps  from 
level  1  through  level  7,  by  the  use  of  equations  (a!7)  and  (al8). 
These  equations  stem  from  an  elementary  re-expressior.  of  the 
horizontal  contour  gradient  at  any  level  as  that  at  a  lovrer  one 
plus  the  gradient  of  the  thickness  between  the  pressure  surfaces 
passing  through  the  two  levels.  The  familiar  hypsometric 
equation  relates  thickness  between  two  pressure  surfaces 
to  the  mean  temperature  in  the  layer.  The  temperature 
gradient  term  on  the  right  band  side  of  (al7),  (al8)  may  be 
viewed  as  a  pseudo-correction  to  be  applied  to  the  horizontal 
contour  gradient  at  a  lower  level  to  obtain  the  horizontal 
contour  gradient  at  the  upper  level.  Through  buildup  from  the 
surface,  only  the  surface  contour  gradient  and  the  horizontal 
temperature  gradients  are  required  as  input. 

It  has  been  assumed  that  the  horizontal  temperature 
gradient  term  is  small  in  the  forest  so  that  up  through  the 
4th  level  the  contour  gradient  is  given  by  the  surface  value 
(equations  alf,  s2£,  j "4  in  si  and  s2).  For  all  levels  above 


31 


4  the  contour  gradient  is  expressed  as  the  surface  contour 
gradient  plus  a  correction  term  which  is  computed  from  equation 
(sl7).  The  value  of  the  surface  contour  gradient  is  generated 
on  the  analog  computer  as  a  linear  interpolation  in  time 
between  the  observed  contour  gradient  at  the  start  of  the 
forecast  period  and  a  predicted  contour  pattern  a;:  the  end  of 
the  forecast  period.  At,  (equation  sl8).  The  Initial  and  final 
contour  gradients  are  scaled  to  100  m/100  Vm  for  contour 
gradients  expressed  in  the  same  units. 

The  horizontal  temperature  gradients  are  routinely 
available  only  at  the  surface  and  at  850  mb.  The  temperature 
gradients  at  the  lowest  and  highest  levels  in  the  forest  and 
the  lowest  and  highest  levels  in  the  free  air  will  be  obtained 
by  linear  interpolaticn  between  the  surface  and  850  mb  levels. 
These  will  be  used  as  input  for  the  solution  and  temperature 
gradients  at  all  other  levels  will  be  obtained  by  linear 
interpolation  between  these.  Thus,  the  correction  term  to 
the  contour  gradient  (equation  al7)  is  expressed  as  a  function 
of  the  horizontal  temperature  gradients  at  the  lowest  and 
highest  levels  (levels  4  and  7  in  our  initial  forest  model). 
Coefficients  and  are  not  simple  linear  geometric 
ratios  but  rather  the  stun  Of  products  of  the  ratio  ~  and 
linear  interpolation  factors.  They  also  contain  conversion 
factors  which  permit  the  temperature  gradients  to  he  expressed 
in  degrees  C/100  km.  The  thickness,  Ah,  is  approximately 
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equal  to  the  height  difference  between  the  two  levels, 
temperature  is  the  mean  for  the  column  between  levels  and 
Zj  but  in  practice  a  constant  temperature  of  298°C  is  used.  The 
percentage  error  involved  in  these  two  approximations  is  very 
small.  The  correction  term  is  scaled  over  a  value  of  A. 34  cm/sec^ 
rather  than  the  3.472  used  for  the  surface  contour  gradient. 
Previous  experience  has  shown  that  the  term  was  underscaled 
if  the  latter  was  used,  for  the  initial  test  solutions  the 
assumption  will  he  made  that  the  horizontal  temperature 
gradient  does  not  change  significantly  over  the  forecast  period. 

2.  Coriolis  Force 

The  coriolis  force  (fu,  fv)  is  carried  over  directly 
into  the  analog  format.  The  computer  variable,  rather  than 
the  physical  variable,  1b  used  in  the  scaled  format.  This  takes 
care  of  the  velocity  scale  factor.  The  time  scale  factor  is 
included  in  the  coefficient  c^,  which  also  incorporates  the 
coriolis  parameter,  f,  a  problem  parameter. 

3.  Kddy-s tress  Terms 

The  eddy-stress  terms  in  the  equations  of  motion  are 
approximated  by  finite  difference  methods,  using  the  values  at 
the  bottom  and  top  of  the  layer  containing  the  particular  level. 
The  assumption  is  made  implicitly  that  the  stress  varies 
linearly  through  the  layer.  The  eddy-stress  terms  are  computed 
at  layer  interfaces,  as  indicated  hv  the  double  subscript 
index  (o.g.  ,  equation  s7,  s8) .  The  nuantity  is  the  x- 
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component  of  the  shearing  stress  at  the  interface  of  the  layers 
containing  levels  5  and  6  (a  height  of  100  ro  in  our  particular 
forest  model).  Similarly  x^  ^  is  the  x-component  of  the  stress 
at  the  interface  of  layers  6  and  7  (200  m  in  our  model) .  The 
coefficient  c„  incorporates  the  density,  the  thickness  of 
the  layer  and  the  scale  factors  for  algebraic  correctness. 

Since  winds  are  not  computed  for  any  levels  above  7,  the 
stress  at  the  top  of  the  layer  containing  level  7  cannot 
be  calculated  directly.  The  assumption  is  made  that  the 
vertical  gradients  of  the  wind  components  at  the  top  of  this 
upper  layer  in  the  free  air  is  the  same  as  it  is  at  the 
bottom.  From  equation  (a7) ,  (aO)  the  ratio  of  the  horizontal 
wind  stress  at  the  top  and  the  bottom  of  the  layer  then 
becones  eoual  to  the  ratio  of  the  momentum  exchange  coefficients 
at  the  top  and  the  bottom  of  the  layer,  that  is, 

Tx,7B  =  Km,78  Ty , 78  =  Km,7S 

Tx ,  67  Nn,67  y ,  67  'm,  67 

This  assumption  permits  x^g  to  be  expressed  in  terms  of  x^ 
in  the  first  equations  of  the  sets  (si)  and  (s2) .  The  factor  . 
which  this  substitution  introduces  is  incorporated  into  the 
coefficient  C£  y 

The  horizontal  wind  stress  for  the  atmospheric  levels  is 
calculated  by  finite  difference  methods  (equations  a7,  a£,  a7f, 
a8f)  assuming  that  the  wind  components  vary  linearly  between 
two  successive  levels,  except  in  the  free  air-canopy  and  forest 
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surface  layers  (a32,  a33,  a30f,  a31f) .  Computation  of  the  stress 
at  the  top  of  the  canopy  is  based  on  the  assumption  that  the  wind 
components  vary  logarithmically  between  the  level  just  below 
and  that  just  above  (a32,  a33) .  The  surface  stress  is  determined 
(a30f,  a31f)  using  an  integral  momentum  exchange  coefficient 
between  the  surface  and  the  first  level. 

The  momentum  exchange  coefficient  at  the  first  level 
above  the  ground  and  at  the  first  level  above  the  top  of  the 
canopy  (eouations  a37f,  n37)  arc  derived  for  wind  profiles 
given  by  Deacon's  8-enuation  and  constant  stress  in  surface 
layer  (0  <  7.  <  A')  and  in  the  upper  air-canopy  layer  (d  < 

7.  <_  d  +  A) .  The  values  of  the  stresses  and  of  the  parameter,  3, 
are  independently  determined  for  the  two  levels.  For  the  initial 
test  simulation  the  exchange  coefficient  is  assumed  to  he 
constant  in  the  free  air-canopy  section  (d  -  A  £  z  5  d  +  A)  and 
to  vary  linearly  between  levels  1  and  3  in  the  forest 
(equation  a!7f).  T.n  the  free  air  above  the  canopy  the  exchange 
coefficient  is  assumed  to  vary  in  accordance  with  the  empirical 
relationship  (al5,  alfi)  developed  during  the  earlier  boundary- 
layer  simulations.  This  assumed  vertical  structure  permits 
the  stress  at  all  other  layer  interfaces  to  be  expressed  in 
terms  of  the  exchange  coefficient  at  the  4th  level  (equations 
s7,  s8,  s32,  s33)  or  as  a  combination  of  the  exchange  coefficients 
at  the  1st  and  4th  levels  (s7f,  s8f)  with  a  considerable 
saving  in  necessary  equipment. 
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The  exchange  coefficients  have  been  scaled  for  50,000  cm2/sec. 
The  coefficient  c^  includes  the  scale  factors  required  to 
maintain  algebraic  accuracy,  the  depth  of  the  layer  separating 
the  two  wind  measurements,  and  the  density,  as  well  as  the 
ratio,  b'^,  between  the  exchange  coefficients  at  the  height  in 
question  and  that  at  the  height  of  level  4.  The  coefficient 
c^^  in  the  stress  equations  for  the  forest  incorporates  the 
depth  between  the  two  wind  measurements  and  the  density. 

Horizontal  Advection  of  Momentum 

The  horizontal  advection  of  momentum  in  the  equations 
of  motion  have  been  re-expressed  in  the  transformation  from 
the  mathematical  to  the  analog  format  so  that  the  horizontal 
gradients  of  the  wind  components  can  be  estimated  from 
kinematic  analyses.  The  advection  terms  become: 

u  (du/dx)  +  v  (du/dy)  •»  a'u  -  8'v 

and  u  (dv/dx)  -4-  u  (dv/dy)  m  a'v  +  0'u 

where  a'  is  the  gradient  of  the  wind  speed  along  the  stream 

line  and  6'  is  the  ratio  of  the  wind  speed  to  the  radius 

of  curvature  of  the  stream  line. 

The  values  of  a'  and  6'  are  routinely  available  only 
at  the  surface  and  the  850  mb  level.  Following  the  method 
outlined  above  for  the  temperature  gradients,  the  values  of 
a'  and  8’  will  be  determined  at  the  lowest  and  highest  levels 
in  the  forest  (levels  1  and  3)  and  in  the  free  air  above  the 
canopy  (levels  4  and  7)  by  linear  interpolation.  These 
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will  be  used  as  input  quantities  for  the  simulation.  The 
products,  n'u,  B'v,  n'v,  ?■  'u  and  the  sum  and  difference 
given  above  are  formed  in  the  simulation  using  the  computer 
generated  values  of  u  and  v  at  levels  1,  3,  4,  and  7,  and  the 
input  a',  O’.  At  all  other  levels,  the  advective  term  is 
obtained  by  linear  interpolation  of  the  advection  calculated 
for  these  four  levels.  Thus  the  advective  term  for  level  6 
appears  as  a  linear  combination  of  the  values  at  levels  4 
and  7  with  the  coefficients,  c and  c^,  representing  simple 
linear  interpolation  factors.  Similarly  the  advective  term 
at:  level  2  appears  as  a  linear  combination  of  the  advective 
terms  at  levels  1  and  3,  In  the  free  air  the  advective 
terms  are  scaled  for  3.472  cm/scc2,  for  the  same  reasons 
that,  were  given  in  the  discussion  of  the  contour  gradient. 

The  scaling  for  the  advective  terms  in  the  forest  (equations 
filf  and  s2f )  is  .6944,  down  by  a  factor  of  five  because  of 
the  lighter  winds . 

5.  The  Drag  Force 

The  analog  transform  for  the  drag  force  in  the  equations 
of  motion  in  the  forest  (slf,  s2f)  is  straightforward. 

The  free  density  and  local  drag  coefficient  have  been  com¬ 
bined  into  a  single  coefficient,  which  is  now  a 

characteristic  of  the  forest,  and  which  will  vary  through 
the  depth  of  the  forest.  The  drag  force  term  has  been 
scaled  for  1000  cm/sec  (equations  slf,  s2f ) ,  but  it  is 
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anticipated  that  this  will  need  to  bf,  changed  following  the 
determination  of  the  magnitude  of  C^.  The  coefficient  c 
is  simply  a  combination  of  scale  parameters  required  to  keep 
the  equations  algebraically  correct. 

Equations  (a9)  and  (alO)  for  determination  of  the 
geostrophic  wind  utilize  the  surface  contour  gradient 
acceleration  term  generated  by  equation  (sl8). 

B.  Air  Temperature  (equations  s3,  s3f) 

The  transformation  of  the  expression  for  temperature 
(equations  3,  3f)  to  the  scaled  analog  format  (s3,  s3f)  is 
similar  to  that  for  the  wind  component  equations. 

1.  Heat  Source  Term,  R 

The  sensible  heat  source  term  which  appears  both  in  the 
free  air  and  in  the  forest  has  been  scaled  over  the  value 
.03472,  the  ratio  of  the  temperature  scale  factor  to  the 
time  scale  factor.  Above  the  canopy  this  term  represents 
primarily  radiational  cooling  arising  from  the  water  vapor 
in  the  atmosphere.  In  the  forest  this  term  i3  more  complicated 
since  not  only  the  atmosphere  but  the  foliage  acts  as  a 
heat  source  (sink)  .  The  manner  in  which  the  forest  term  will 
be  handled  has  mt  been  settled  and  is  under  investigation. 
Above  the  canopy  the  approach  will  be  similar  to  that 
which  has  been  employed  in  previous  simulations  on  the  GPAC. 

2.  Convective  Heat  Flux 

The  convective  heat  flux  terms  are  analogous  to  the  stress 
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terms  in  the  momentum  equations  and  analog  forms  for  these 
terms  in  the  temperature  equations  (a3  and  a3£)  are  very 
similar  to  the  analog  forms  of  the  stress  terms.  Since  the 
temperature  is  not  computed  for  levels  above  7,  the  flux  at 
the  upper  boundary  of  the  top  layer  is  expressed  as  a  function 
of  that  at  the  lower  boundary  by  using  the  same  argument 
that  was  used  for  the  horizontal  stress  at  the  same  boundary. 

That  is,  the  vertical  temperature  gradient  is  assumed  to  be  the 
same  at  the  bottom  and  top  of  the  uppermost  layer  in  the 
free  air  section  so  that,  from  equation  (a9),  the  ratio  of  the 
flux  at  the  top  of  the  layer  to  the  flux  at  the  bottom  of  the 
layer  is  equal  to  the  ratio  of  the  heat  exchange  coefficients 
at  the  top  and  bottom  of  the  layer.  Thus  the  calculation  of 
the  temperature  at  level  7  (equation  s9)  involves  only  the 
convective  flux  at  the  interface  of  levels  6  and  7. 

The  convective  heat  fluxes  are  obtained  from  equations 
(s9)  and  (sllf),  under  the  implicit  assumption  that  the 
temperature  varies  linearly  between  two  successive  levels. 

Since  the  heat  exchange  coefficient,  K^,  is  assumed  equal 
to  the  momentum  exchange  coefficient,  the  value  at  any  level 
above  the  forest  can  be  expressed  as  a  function  of  the 
exchange  coefficient  at  level  4  by  equations  (al5)  and  (al6). 

In  the  forest  the  exchange  coefficient  is  given  as  a  function 
of  those  at  the  first  and  fourth  levels  by  using  the  relationship 
given  in  equation  (37). 
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Since  the  vertical  temperature  gradient  just  above  the 
canopy  may  be  different  from  that  in  the  upper  portion  of  the 
canopy  the  upward  convective  flux  at  the  air-canopv  interface 
is  calculated  separately  from  the  downward  flux  from  the  same 
level.  The  former  is  based  on  the  change  of  temperature 
in  the  small  layer  just  above  the  interface  (equation 
s25)  and  the  latter  on  the  change  of  temperature  in  the  small 
layer  just  below  (equation  s26). 

3.  Horizontal  Temperature  Advection 

As  might  be  expected  the  temperature  advection  terms 
are  handled. in  exactly  the  same  manner  as  the  wind  advection 
terms.  That  is,  the  temperature  gradients,  3T/9x  and  3T/3y, 
like  a  and  6,  are  obtained  at  levels  1,  3,  4,  and  i .  by 
linear  interpolation  between  surface  and  850  mb  and  only 
these  are  used  as  input.  From  this  input,  [u(3T/Dx)  and 
v(3T/8y) ]  at  the  same  levels  are  computed  and  linear  interpolation 
is  employed  to  obtain  the  values  at  the  other  levels. 

The  values  of  the  horizontal  temperature  gradients  must 
be  expressed  in  units  of  deg/ cm  for  this  equation  to  be 
used  directly  since  the  coefficients  c^and  are  simply 
interpolation  factors.  For  the  initial  simulation  the 
horizontal  teippe^ture  gradient  will  be  treated  as  invariant 
in  time. 

C.  Vapor  Pressure 

The  defining  equations  for  the  atmospheric  vapor  pressure 
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(equations  4  and  4f)  are  very  similar  to  those  for  air  temperature 
(3  and  3f ) .  Not  too  surprisingly  th  analog  and  scaled  forms 
are  also  very  similar  to  those  for  temperature  and  were 
obtained  in  a  manner  almost  identical  to  that  described  in 
section  II  B,  above. 

The  vapor  source  term,  M,  in  the  forest  (s3f)  includes 
both  evapo- transpiration  of  the  vegetation  and  condensation 
(evaporation)  in  the  atmosphere.  For  the  initial  simulation  the 
latter  will  be  assumed  to  be  zero  in  the  free  air  above  the 
forest,  i.e.,  none  will  be  permitted. 

In  the  transformation  from  mathematical  (equations  10  and 
12f)  to  analog  format,  the  evaporative  heat  flux  has  been 
re-expressed  in  terms  of  the  vapor  pressure  utilizing  the 
relationship  given  in  equation  (11).  This  substitution 
introduced  the  pressure  at  each  simulation  level  which,  under 
the  hydrostatic  assumption  and  the  assumption  of  constant  density, 
are  determined  from  equation  (a5) .  As  in  the  case  of  the  stresses 
and  convective  fluxes  the  solution  of  the  evaporative  fluxes 
is  based  on  finite  differencing  techniques.  In  addition  there  has 
been  some  algebraic  manipulation  in  arriving  at  equations  (slO) 
and  (sl2f)  so  that  only  one  of  the  factors  in  the  difference  on 
the  right  hand  side  contains  a  multiplier  and  the  pressure  at 
only  one  level  appears  in  the  equation  explicitly.  Beyond 
this  the  evaporative  heat  fluxes  are  obtained  in  a  manner  similar 
to  that  described  for  the  convective  fluxes  not  ’’nl.y  in  the 
free  air  and  forest  sections  but  also  at  the  air-canopy  interface. 
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III.  Diffusion  experiments  in  the  Panamanian  Forest 

Smoke  releases  were  made  in  the  forest  of  the  Panama  Canal 
Zone  during  the  last  week  of  August.  The  chief  purpose  of  these 
releases  was  to  familiarize  contract  personnel  with  the  air  flow 
in  a  tropical  forest.  Thus  the  tests  were  qualitative  in  nature 
and  observations  were  primarily  visual  or  photographic.  A 
concerted  effort  was  made  to  obtain  some  idea  as  to  the  influence 
of  terrain,  tree  density,  time  of  day,  cloudiness,  local 
clearings,  and  surface  brush. 

Several  sites  having  different  terrain,  tree  types  and 
foliage  characteristics  were  selected  for  the  tests.  Although 
the  density  of  the  forest  at  the  sites  varied  considerably, 
the  depth  of  the  forest,  i.e.,  the  height  of  the  top  of  the 
canopy,  was  fairly  uniform  around  100  to  l?0  ft.  Smoke  pots 
and  grenades  provided  the  tracer  material.  The  grenades  pro¬ 
vided  smoke  releases  in  ravines,  streams,  slopes,  and  generally 
inaccessible  locations.  They  were  nearly  always  used  in 
conjunction  with  smoke  pots  located  at  more  level  locations. 

All  release  points  were  on  the  ground. 

Very  little  supporting  meteorological  data  are  available. 
Three  tests,  at  one  hour  intervals  on  a  single  day,  were  made 
at  the  Albrook  Data  Base  Site.  For  these  tests,  recorded  measure¬ 
ments  are  available  at  5  minute  intervals  for  the  wind  direction 
and  speed  at  the  46  m  level  (about  16  m  above  the  forest  canopy) 


and  for  temperature  and  dewpoint  at  a  number  of  levels  on  the 
tower.  The  wind  measurements  in  particular  are  valuable  since 
they  provide  information  about  the  amount  of  turning  that  occurs 
through  the  canopy.  For  all  but  two  of  the  other  tests, 
tethered  pilot  balloons  were  used  to  estimate  the  direction  of 
the  wind  above  the  canopy.  This  method  was  not  very  satisfactory. 
Our  limited  observations  are  briefly  summarized  below. 

1.  There  was  unmistakable,  and  frequently  pronounced,  horizontal 
transport  of  the  materials  in  the  forest  and,  except  for  an 
initial  vertical  motion,  this  was  the  dominant  transport. 

2.  There  was  spreading  of  the  plume,  indicating  diffusion  of 
material  around  the  mean  direction  of  the  flow.  The  behavior 
of  the  smoke  was,  generally,  not  very  different  from  that  of 
stack  plumes  with  moderate  flue  temperature  at  more  open  sites. 

3.  The  horizontal  transport  appeared  to  be  nearly  constant  in 
direction  from  about  four  meters  above  ground  to  at  least  sixty 
to  seventy  percent  of  the  height  of  the  forest, 

4.  At  Albrook  the  direction  in  which  the  smoke  plume  moved 
through  the  forest  was  within  30°  of  the  wind  direction  16  m 
above  the  forest.  It  was  always  to  the  left  of  the  wind, 
suggesting  clockwise  turning  (veering)  of  the  wind  with  height 
through  the  canopy.  In  many  of  the  other  tests  also,  the 
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transport  in  the  forest  appeared  to  be  20  to  45  degrees  to  the 
left  of  the  wind  above  the  forest.  However,  the  tethered 
balloon  estimates  of  the  free  air  wind  are  not  very  accurate. 

5.  Although  the  horizontal  transport  tended  to  maintain  a  mean 
direction,  oscillations  in  direction,  with  periods  of  a  few 
minutes,  were  noticed  in  nearly  every  test.  The  oscillation 
usually  had  an  amplitude  of  about  ±  30°  of  the  mean  transport 
direction,  but  at  least  once  it  was  more  pronounced  (directional 
change  of  130-150  degrees)  and  also  of  longer  period  (eight  to 
ten  minutes).  The  Albrook  observations  16  m  above  the  top  of  the 
canopy  indicated  variations  in  wind  direction  of  the  same  magnitude 
as  that  observed  for  the  smoke  plume  in  the  forest.  Because  of 
the  relatively  long  time  between  observations  (5  minutes) ,  it 

is  difficult  to  ascertain  the  period  of  the  directional  deviation. 

6.  At  only  one  site  was  smoke  clearly  seen  above  the  canopy. 

This  occurred  in  the  thinnest  of  the  forest  sites  in  early 
afternoon.  There  may,  however,  have  been  some  movement  up 
through  the  canopy  in  other  tests  where  poor  viewing  conditions 
made  it  difficult  to  detect  the  smoke  in  the  free  air. 

7.  Diurnal  changes  could  be  checked  on  only  the  one  day  when 
three  runs  were  made  at  hourly  intervals  between  0930  and  1130, 

The  horizontal  transport  was  predominantly  from  the  south¬ 
west  in  the  first  test,  from  the  east  on  the  second  and 
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principally  from  the  west  in  the  third.  These  tests  were  made 
at  the  Albrook  Tower  Site  and  the  recorded  winds  16  m  above  the 
forest  show  the  same  longer  period  variation. 

8.  Terrain  appears  to  affect  the  horizontal  transport  very 
little.  The  smoke  from  grenades  thrown  into  stream  beds  and 
gullies  tended  to  move  up  and  out  and  take  on  the  motion  of  the 
general  horizontal  transport  as  indicated  by  the  smoke  from 
releases  at  more  level  sites  in  the  vicinity.  Similarly  smoke 
grenades  thrown  up  on  to  slopes  tended  to  drift  in  the  general 
direction  of  the  horizontal  transport  and  did  not  tend  to  flow 
either  up  slope  or  down  slope  unless  they  coincided  with  the  mean 
direction.  Small  hills  or  knolls,  however,  did  appear  to  cause 

a  decrease  in  wind  speed  in  the  lowest  5  m  immediately  downwind. 

9.  Local  clearings  and  roads  seemed  to  have  little  effect  in 

the  general  transport  of  the  smoke.  The  plume  from  forest  releases 
were  observed  to  cross  clearings  and  wide  roads  and  to  re-enter 
the  forest  on  the  other  side.  However,  the  cross-wind  diffusion 
of  material  was  enhanced  in  the  clearings  and  particularly  at 


wide  roads. 
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APPENDIX  A 
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APPENDIX  B 


A.  Free  Air  Section 
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